Pore dynamics of osmotically stressed vesicles 
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We present a theory for pore dynamics of osmotically stressed vesicles. When a liposome with 
an internal concentration of solute is placed inside a solute-depleted medium, an osmotic flow of 
solvent through the lipid bilayer leads to swelling of vesicle and to increase in membrane surface 
tension. This can result in membrane rupture and opening of thermal pores. Depending on the 
internal concentration of solute and the size of the vesicle, pores can close rapidly or be long-lived. 
We find that the life span of the long-lived pores scales non-trivially with the size of the liposome. 
Closure of the long-lived pore is followed by a rapid flicker-like opening and closing of short-lived 
pores. Our model is consistent with the observation of long-lived pores in red blood cell ghosts. 



Opening of stable pores in cellular phospholipid mem- 
branes is an important step for drug delivery and gene 
therapy. A number of methods have been explored to this 
end. Electroporation [jj is the method of choice when 
gene delivery is performed in vitro. It is, however, associ- 
ated with a significant cell mortality limiting its practical 
use. Holes in cellular membranes have also been opened 
by intense illumination 0, The mechanism leading 
to pore formation under these circumstances is still not 
fully understood. It, however, appears to be similar to 
the action of laser tweezers Q , a part of the low dielec- 
tric constant lipid membrane is sucked into the laser trap, 
stretching the membrane until rupture. 

For drug delivery the goal is to introduce a controlled 
amount of substance to the specific disease site. The de- 
livery system should minimize the potential side-effects 
by selectively targeting the affected region of the organ- 
ism [j| . Liposomes can be made to contain drugs or genes 
and can be introduced into blood stream without pro- 
voking an immune system response. Liposomal affinity 
for specific tissue can be increased by varying the mem- 
brane lipidic composition or by including ligands that 
recognize and bind to specific receptors |6|]. One exam- 
ple of this is the drug called Doxirubicin used in cancer 
chemotherapy. The liposomal encapsulation of Doxiru- 
bicin has been shown to efficiently deliver it to tumors, 
while minimizing the side effects, such as damage to heart 
muscle tissue. 

When liposome ruptures, its content leaks out. It is 
precisely the dynamics of this rupture that is the subject 
of the present Letter. In particular, our goal is to un- 
derstand the mechanism of pore opening and closure in 
the osmotically stressed vesicles. The situation which we 
attempt to model is the following. Liposomes are pre- 
pared inside a solution containing small molecules im- 
permeable to lipid membrane at concentration cq. They 
are then placed into a purified aqueous solution. The 
imbalance in solute concentration inside and outside the 
vesicles leads to osmotic flow of solvent through the semi- 
permeable lipid membrane, resulting in swelling of vesi- 
cles and stretching of their membranes. As the ten- 
sion increases, the energy barrier to pore nucleation de- 
creases |3- If there is enough solute inside the lipo- 
some, the barrier height eventually becomes comparable 



to the energy of thermal fluctuations. When this hap- 
pens, membrane raptures and a pore opens One in- 
teresting question is what is the life span of this pore? 
Specifically, as a pore opens it relaxes the membrane ten- 
sion. However, opening of a pore exposes the hydropho- 
bic interior of membrane to water, leading to pore line 
(edge) tension. This hydrophobic effect induces a force 
which tries to re-seal the pore. Furthermore, a hole in 
membrane results in efflux of both solute and solvent, re- 
ducing the vesicle volume and surface tension. All these 
effects lead to an intricate interplay of time and length 
scales which control the pore closure. 

We shall designate the radius of a vesicle as R and the 
radius of a pore as r. The conservation of mass then 
leads to 

4npR 2 ^ = j w - Trr 2 pv , (1) 



dt 



where p is the density of water, j w is the osmotic flux and 
v is the leak-out velocity. The leak-out velocity 0, 0, 
ITU IT^| is determined by the balance between the shear 
stress, proportional to r/v/r, and the osmotic pressure 
inside the vesicle Ap, 
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where rj is the solvent viscosity. Inside a swollen vesi- 
cle the osmotic pressure is compensated by the Laplace 
pressure, so that 

Ap=|, (3) 

where a is the tension of a stretched membrane. In this 
work we shall concentrate on vesicles with fairly large 
internal concentrations of solute, close to 1 M. This cor- 
responds to osmotic pressures as high as 20 atmospheres. 
Of course, before these extreme pressures can be built, 
the lipid membrane will rapture, leaking out some of the 
internal content and releasing the surface stress. At these 
high pressures, the thermal undulations of membrane can 
be ignored, and the membrane tension is controlled by 
the stretching modulus [Tiil H3. Hl| K s . The membrane 
elastic energy takes a Hook-like form 

E. = l 
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where Aq = AttRq is the equilibrium surface area of the 
unstretched membrane. As the membrane is stretched 
beyond the elastic limit, pores are nucleated reducing 
the membrane tension. The elastic energy of a membrane 
with a pore is 
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E{R,r) = -K s [An(R 2 - R 2 
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where r is the pore radius and 7 is the pore line tension. 
The membrane surface tension is, 



/ p s dE s 



K s (Att(R 2 -R 2 )~nr 2 ) . (6) 



The osmotic current j w is determined by the membrane 
permeability P, the concentration difference of solute in- 
side and outside the vesicle, and the Laplace pressure. A 
simple phenomenological expression for the osmotic cur- 
rent of water into the vesicle can be written as, 
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where ks is the Boltzmann constant and Na is the Avo- 
gadro number. If the difference of molar (M) solute 
concentration c, inside and outside the vesicle is not 
very large, the integrity of membrane will not be com- 
promised, and a stationary state with j w = will be 
achieved. Under these conditions the osmotic pressure 
is completely compensated by the Laplace pressure, re- 
sulting in zero net flux of solvent. For sufficiently large 
internal concentration cq, a stationary state will not be 
achieved before membrane ruptures. It is precisely in this 
regime that we expect to see some interesting physics. 

The growth of a pore is controlled by the rate at which 
the elastic energy is dissipated by the membrane viscosity 
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where I is the membrane width. Since the membrane is 
impermeable to solute particles the internal solute con- 
centration is modified only through the osmotic influx of 
solvent and the efflux of solute through the open pore, 
after membrane has ruptured. The continuity equation 
expressing this is 
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where we have assumed that solute is uniformly dis- 
tributed inside the vesicle. 

We are now in position to study the evolutionary dy- 
namics of an osmotically stressed vesicle. The char- 
acteristic values of the parameters appearing in the 
above equations are: AttRqK s = 0.2 J/m 2 , I — 3.5nm, 
7 = I CP 12 J/m, r) m = 100 Pas , rj w = 10~ 3 Pas and 
P = 1.8 X W- i kg/m 2 sM. 

The energy necessary to open a pore or radius r is 
AE{R,r) = E(R,r) - E(R,0). The form of AE(R,r) 
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FIG. 1: Energy AE(R,r) necessary to open a pore of radius 
r in a liposome with Ro = 100 nm. The curves are, from 
top down, R/R = 1.0, R/R = 1.0008, R/R = 1.001, and 
R/Ro = 1.0012 



as a function of pore size is plotted in Fig. 1 for vari- 
ous values of R. As the osmotic flow of solvent stretches 
the membrane, a minimum in AE(R, r) for r ^ indi- 
cates energetic favorability of pore formation. However, 
to reach this minimum a sufficiently large thermal fluc- 
tuation is necessary to overcome the energy barrier. The 
probability of such a fluctuation occurring is proportional 
to the Boltzmann factor 



P(r) 



(10) 



where (3 = 1/ksT and r* is the critical pore size. Pores 
with r < r* will shrink and close while pores with r > r* 
will grow. The waiting time for appearance of a thermal 
fluctuation with enough energy to open a pore of radius 
r* is very long, unless AE(R,r*) w fc^T. 

In the limit of large surface tensions the critical pore 
size is r* w 7/cr and the barrier height is 
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We note that when membrane tension reaches a c = 
10 -3 J/m 2 , the barrier height is approximately fc^T. 
This value of a c is in agreement with the tensions found 
to be necessary to rupture a mechanically stretched mem- 
brane p"?l Il8j| . The minimum solute concentration nec- 
essary for a pore to open is 
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which for a vesicle of 100 nm is approximately IQmM. 

The dynamics then proceeds as follows: Eqs. (0,10, 
and are solved numerically, using Euler's method with 
a stepsize of 8t = 0.1/is, to find the evolution of R(t), 
r(t), and c(t). When the surface tension becomes such 
that, AEb w ksT, a pore of size r* opens, see Fig. 2, 
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and the internal content of the liposome begins to leak- 
out. This process reduces the volume of the vesicle and 
decreases the membrane tension, until the pore closes. 
Once this happens, the osmotic swelling of the liposome 
re-starts. As the energy barrier to pore nucleation drops 
down to ksT a new pore is opened, etc. The process 
stops when the internal solute concentration reaches c™ m 
and a stationary state with j w = is established. 
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FIG. 2: Radius of a pore as a function of time for vesicles of 
Ro = 200 nm (topmost), 220 nm and 240 nm (bottommost), 
for initial concentration Co = 0.5M. Note that for vesicle 
with Ro — 200 nm the pores are short-lived, while for larger 
vesicles, long-lived pore opens first. 

For small vesicles, pores open and close very quickly 
with the characteristic time r/ ~ 1CP 2 , resulting in a pe- 
riodic flickering. For liposomes with Rq > -Rg(co) a long- 
lived pore appears, see Fig. 3. The critical size of a lipo- 
some Rq(cq), necessary for nucleating a long-lived pore 
depends on the internal solute concentration. Larger the 
initial solute concentration cq, smaller will be the size 
of the vesicle which supports a long-lived pore. The long 
life span of these pores is the result of a "wash-out" effect 
in which the osmotic flux is almost completely compen- 
sated by the leak-out rate of solute through the pore. 
When solute concentration inside the vesicle drops below 
the critical value, c c (Rq), the long-lived pore closes. This 
value is insensitive to the initial solute concentration cq, 
but depends strongly on the vesicle size Ro. The life span 
of a long-lived pore r scales with the vesicle size Rq, 



(13) 



with v sa 2.3 — 2.4, see Fig. 3. After the long-lived pore 
has closed, it is followed by a sequence of short-lived pores 
with the characteristic life span 77. 

We have checked that the pore dynamics is not sensi- 
tive to the specific mechanism of pore opening. In par- 
ticular, even if the pores are opened stochastically, with 
the probability given by the Boltzmann distribution, the 
rate of leakage and the life span of the long-lived pores 
is affected very little. 
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FIG. 3: Life span of the first-open pore as a function as a 
function of vesicle size Rq for Co = 0.5M. Note the appearance 
of the critical vesicle size Rq(co) which can sustain long-lived 
pores. 



Up to now we have not taken into account the diffusive 
efflux of solute through the open pores. The diffusive 
current through a pore of size r is approximately 



3d ~ 7rr c 
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Ro 
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This leads to the decline of internal solute concentration 
governed by the conservation equation, 



47r 3 dc 

the solution of which is 

c(t) = c Q e~ 



-3d 



where the effusion time r e is 
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3r 2 D 
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Using D w 10 _9 m 2 /s, appropriate for small organic 
molecules such as sucrose, and r = r* ss j/a c « 1 nm, 
we see that for liposomes with R = 50 nm, the time 
for effusion is r e « 1CU 2 s. This is comparable to the 
flicker time 17. Therefore, for small vesicles effusion is 
an important mechanism for the loss of solute. On the 
other hand, for liposomes with Rq = 100 nm and above, 
effusion is only marginally relevant. 

Lon g lived p ores have been observed in red blood cell 
ghosts 19, 20], their size dependent on the ionic strength 
of the surrounding medium |21| . No theory, up to date, 
was able to account for these long-lived pores. Holes 
were predicted to either grow indefinitely, which would 
result in ghost vesiculation, or to close completely |22^ . 
Our model provides a dynamical mechanism for pore sta- 
bilization, consistent with the experimental observations. 
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In the specific case of red blood cell ghosts the ratio j/K s 
has to be ajusted to account for the large radius r/R ob- 
served in these experiments 0, ■ 

In aqueous solutions the phospholipid membranes ac- 
quire a net negative charge. At physiological concentra- 
tions, 154 mM of NaCl, the Debye length, however, is 
quite short, less then 1 nm and the electrostatic inter- 
actions are strongly screened |23|. We, therefore, do not 
expect that electrostatics will significantly modify the ba- 
sic conclusions of our theory, beyond the renormalization 
of membrane line and surface tension. However, fur- 
ther, investigations in this direction are necessary and 
will be the subject of future work. 

Finally, it is curious to note a strong similarity be- 
tween the rupture of osmotically stressed liposomes and 
bursting of Hydra cells aggregates |24J . Hydra and marine 
sponges can generate functional organisms from random 
cell aggregates purely through the intercellular interac- 
tion. This morphogenesis is characterized by the cav- 
ity formation followed by swelling and violent bursting, 
which expels the internal fluid and dead cells, see Fig. 5. 



FIG. 4: A sequence of bursts in hydra aggregates. Picture one 
shows an aggregate with a diffuse cloud of expelled cells from 
the previous explosions, while picture two shows a cellular 
aggregate in the process of bursting. Arrow indicates the site 
of the burst. 



Unlike the liposomes, cellular aggregates are too mas- 
sive for temperature to be of any relevance. Instead the 
cellular Brownian motion is driven by the metabolic fluc- 
tuations of the cytoskeleton j^. We can, therefore, ex- 
pect that a theory similar to the one presented above 
for the liposomes, might also apply to osmotically driven 
bursting of cellular aggregates. 
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